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In this talk we would like to review recent results on non-singular cosmological models. It has 
been recently shown that among stiff perfect fluid inhomogeneous spacetimes the absence 
of singularities is more common than it was expected in the literature. We would like to 
generalize these results and apply them to other matter sources. 



1 Introduction 

The interest for non-singular cosmologies has risen since it was shown that among inhomo- 
geneous stiff fluid spacetimes there was a regular family depending on two nearly arbitrary 
functions [1]. 

Since the publication of the first known non-singular perfect fluid cosmological model by 
Senovilla [2], the number of regular spacetimes had not increased a lot. Senovilla and Ruiz 
enlarged the family [3] two years afterwards. All of them had the equation of state of an 
incoherent radiation fluid. The other regular models that were found either had no equation 
of state [4] or were for stiff fluids [5, 6, 7, 8]. 

As it has been said, the role of the stiff fluid equation of state in the absence of singu- 
larities is well known, since the gradient of pressure prevents the formation of singularities. 
In fact, it is a limiting case between causal and non-causal linear equations of state, since 
the velocity of sound for the stiff fluid is equal to the velocity of light. On the contrary, dust 
models, with vanishing pressure, are all singular. 

Therefore, Senovilla's original solution is intriguing for the fact that it contains incoherent 
radiation as a matter content, whereas every other barotropic regular model has been found 
for a stiff equation of state. Furthermore, in the generalization of the model by Senovilla 
and Ruiz incoherent radiation spacetimes are the only ones that are non-singular and they 
trace the separatrix between models which have Ricci and Weyl singularities. 

In this talk we pursue to analyze the role of radiation in the absence of singularities 
in inhomogeneous cosmological models. To this aim we shall consider pure radiation cos- 
mological models in the framework of orthogonally-transitive, mutually orthogonal Abelian 
G2 spacetimes, where most regular perfect fluid models have been found so far. There are 
spherically symmetric non-singular models [9] , but the matter content is an anisotropic fluid. 

In the next section we shall derive the equations for pure radiation inhomogeneous cos- 
mological models and in the last section we shall analyze the geodesic completeness of such 
models. The results will be summarized in a conclusions section. 
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2 Pure radiation inhomogeneous models 

We start by defining the matter content of the models. The energy-momentum tensor for 
pure radiation, 

Ti"" = ^kf^k", (1) 

is written in terms of a null vector field describing the 4- velocity of the radiation field k and 
a positive function ^. It is tempting to identify with the energy density of the radiation 
field, but it must be noticed that it is not unambiguously defined, since it may be altered 
by rescaling the vector field k. 

In full generality, conservation laws for this energy-momentum, 

T';^ = ^^^k^k" + ^a" + ^Oki" = 0, (2) 

are written in terms of the expansion = kY^ and the acceleration, = ki\,k" of the vector 
field k. 

If k were timelike instead of null, we could split the equation in two components, parallel 
and orthogonal to k. But in the null case there is a single component if is non-zero. 

As we have already stated, we require that the spacetimes allow an orthogonally- 
transitive, mutually orthogonal G2 group of isometrics. The metric, 

ds^ = e^^'i-dt^ + dx^) + p2e2C^dy2 + e-^^'dz^, (3) 

is written in coordinates y, z, adapted to the isometry generators, which do not appear in 
the metric functions p, K, U, and two isotropic coordinates t and x. The function K is just 
a conformal factor and p is the surface element of the transitivity surfaces. It is much related 
to the absence of singularities, since all regular models have a spacelike gradient of p. The 
coordinates have the usual ranges. 

We impose that k has just components on t and x, 

k = e-''{dt + d,), (4) 

for an outgoing radiation flux. For ingoing radiation we have just to change the plus sign 
for a minus in the sum. 

In these coordinates the conservation equation can be integrated, 

{dt + d,) {^e^^p} = ^ ^e^^p = h{t - x), (5) 

in terms of one function h of the outgoing null coordinate t — x. 
The remaining Einstein equations in these coordinates, 

Utt-Ua;x + -{Utpt-UxPx)=0, (6a) 
P 

Ptt - Pxx = 0, (6b) 

KtPx + K^pt = ptx + UtPx + UxPt + 2pUtUx - h, (6c) 

KtPt + K^Px = + Utpt + U^Px + P {Ui + C/J) + h, (6d) 

may be reduced to a plane-wave equation for U and a wave equation for p. The equations 
for K are a quadrature. The intcgrability condition for these equations is a consequence of 
the other equations and therefore they can be integrated at the end. 

If we take p as a coordinate, the system is further simplified. There are three possibilities: 
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If p is a null coordinate, p = p{t — x), U = U{t — x) is also null and the only equation is 
a wave equation, 

Ktt - i^.. = ^ K{u, v) = f{u) + g{v), (7) 

which may integrated in terms of null coordinates, u = t — x, v = t + x. Just / is 
constrained by the quadrature, 



- — \-U + 1 

2p' p' 2p' 

where the prima denotes derivation with respect to u. 

After eliminating redundant functions by defining u = J due^^ , v = J dve^^, the metric 
is written in a simple form, 

ds'^ = -2dudv + F'^dy^ + G^dz^, (9) 

in terms of two arbitrary functions, F and G, of u. We have dropped the tilde for 
simplicity. 

The function ^ must be positive, 

These models have a five-dimensional group of isometrics acting on null orbits. The vector 
field dv is covariantly constant and these spacetimes contain plane-fronted gravitational 
waves with parallel rays. 

• If p is a spacelike coordinate, we choose p = x, 

Utt -U^.- — =0, (11a) 

X 

Kt = Ut + 2xUtU^ - h, (lib) 
K, = U^+x{U^ + Ul) +h, (11c) 

and the equations are the same as for vacuum spacetimes except for the function h, which 
appears in the metric as a conformal factor —H{t — x) added to K [10]. 
Since according to equation 5 the fimction h{t — x) must be positive for positive x and 
negative for negative x we have to interpret a; as a radial coordinate ranging from zero 
to infinity, although the axis is not flat. Therefore Ht = h must be positive and H is a 
growing function. 

• Finally, if p is a timelike coordinate, we take p = t, 

Utt - C/x. + y = 0, (12a) 

K., = U^ + 2tUtU^-h, (12b) 
Kt = Ut+t (f/2 + C/2) + /i, (12c) 

and the results are quite similar to the spacelike case. We shall not develop this case 
since it is singular at f = 0. 



3 Geodesic completeness 



There are at least three different definitions for singularities [11]. One of them states that a 
spacetime is singular if any of the polynomial curvature scalars blow up at an event. This 
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definition has the drawback that singular events may be hidden by conveniently removing 
an open set aroung them. 

Therefore a more reasonable definition would be considering causal observers on tra- 
jectories on the spacctime. For instance, restricting to observers in free fall, which move 
along timelike and null geodesies, a singularity is encountered if their proper time cannot 
be extended from minus infinity to infinity. This would mean that the observer has left the 
spacetime in a finite proper time. Spacetimes which do not suffer this inconvenience are 
called causally geodesically complete. We shall adopt this definition for regular spacetimes. 

There is another definition of singularity which refers to general causal observers, not 
necessary in free fall and therefore not following causal geodesies. 

In order no analyze the regularity of our spacetimes we are lead to study their causal 
geodesies. 

A geodesic of 4- velocity u, 

u^ = i^=^^, (13) 
where r is the proper time, satisfies the geodesic equations written in the form, 

a^^ = x^ + ri^pX^xP = 0, -00 < r < 00, (14) 

in terms of the Christoffel symbols for the metric of the spacetime. 

Since these models have isometrics, the geodesic equations become simpler, because there 
are conserved quantities of geodesic motion assigned to them. If ^ is a generator of an 
isometry, the quantity X, 

X := g^.eW, (15) 
is conserved along geodesies of 4- velocity u, 

DuX = X^^u'' = (.g^,.;.r«," + + 9>^,^C^^-a) = 0, (16) 

since the metric g is covariantly constant, the acceleration a'^ = u^^u" of the geodesies is 
zero and the generator satisfies the Killing equation, 

W+^-;^=0. (17) 
Besides, the character of the geodesic is preserved. 



(5 -S^.M^u", (18) 

and equal to zero for mdl geodesies and to one for timelike geodesies. In fact, geodesic 
equations merely state this fact. 

All these conserved quantities will be useful. 

We have classified the spacetimes by the character of the coordinate function p: 

• The function p is a null coordinate: We have seen that the vector fields {9^,, 9y, are 
generators of isometrics and therefore provide constants of geodesic motion. 



Y = {dy,u)=g^^u^ = F-^V, 
Z={d,,u)=g,^u^ = G''z. 



(19) 
(20) 
(21) 
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All these equations are integrable provided that F^^ and are integrable functions. 
We are left just with one equation for v, which we may derive directly from the expression 
for 5, without using the Christoffel symbols, 



These results are comprised in a simple theorem: 

Theorem 1. A spacetime endowed with a m,etric (9) and with a pure radiation field, as 
matter content is causally geodesically complete if and only if the metric components g^^ 
and g^^ are integrable functions of the null coordinate u. 

However, as it is explained in [12], the positiveness requirement for the function $ implies 
that the metric functions are continuous hut not differentiablc. 
• The function p is a spacelike coordinate: In this case there are in general no additional 
isometries. However, there are theorems specifically derived for this case [13]; 

Theorem 2. A diagonal Abelian orthogonally transitive spacetime with spacelike orbits 
endowed with a metric in the form (3) with metric functions K,U,p, where p has 
a spacelike gradient, is future causally geodesically complete provided that along causal 

geodesies: 

1. For large values of t and increasing x, 



a) {K - U - \np)t + {K - U - Inp)^ > 0, and either {K - U - Inp)^ > or 

\{K-U- \np),\ < (K -U-\np)t + {K-U- Inp),. 

b) Kt + K^> 0, and either K^>0 or \K^\ <Kt + K^. 

c) {K + U)t + {K + U)^ > 0, and either {K + U)^>Q or\{K + U)^\ < (K + U)t + 



For geodesies pointing to the past the theorem imposes the same conditions but reversing 
the sign in time derivatives. 

This result imposes certain restrictions on the metric functions: 

Theorem 3. A diagonal Abelian orthogonally-transitive spacetime with spacelike orbits 
with a metric in the form (3) and a spacelike p and pure outgoing radiation as matter 
content is causally geodesically complete if: 

1. x^~^\Ux ± ?7i| 7^ for large values of \t\ and x. 

2. H{t) < In \t\ + b and U(t, 0) > 6 - i In jtj + H{t) for large values of t. 

3. U{t, Q)>b-\\n\t\ + H{t)/2 for small values oft. 

This result has obvious consequences. The amount of energy that can be radiated, which 
may be measured by H , must be bounded in order to prevent the formation of singular- 
ities. For ingoing radiation the result is very similar. 

Theorem 4. A diagonal Abelian orthogonally transitive spacetime with spacelike orbits 
with a m,etric in the form, (3) and a spacelike p and pure ingoing radiation as matter 
content is ca,usa,lly geodesically complete if 

1. x^~'^\Ux ± t^tl 7^ for large values of \t\ and x. 

2. U{t, 0) > 6 - i In jtj - H{t)/2 for large values oft. 

3. H{t) > - In jtj + fe a,nd U(t, 0) > 6 - i In |t| - H{t) for small values of t. 

Therefore we may derive regular radiation cosmological models with a function U which 
grows for large values oi \t\ and x and a function H that does not grow or decrease too 
quickly. 



dv _ 1 
d^~2V^ 



{5 + Y^F-^ + Z^G-^}. 



(22) 




6 L. Fernandez- Jambrina and L.M. Gonzalez-Romero 

4 Conclusions 

We have derived sufficient conditions for an inhomogeneous spacetime filled with pure radi- 
ation to be singularity-free in the sense of causal geodesic completeness. The conditions are 
fairly easy to implement in terms of a growing function H and a solution to the plane- wave 
equation with very few restrictions. This is much the same as it happened for stiff perfect 
fluids [14]. Regular models are pretty general among radiation spacetimes. 

Pure radiation fields with positive energy density fulfill weak, strong and dominant energy 
conditions. Since these models have a function t with timelike gradient everywhere, they are 
causally stable [11] and hence they satisfy every weaker causality condition. For instance, 
they do not contain closed causal curves. The only possibility to circumvent the powerful 
singularity theorems is the absence of closed trapped surfaces. 

It has already been shown that singularity-free models are common among stiff perfect 
fiuids and pure radiation inhomogeneous spacetimes. It remains to be seen if these results 
can be extended to other matter contents and symmetries. 
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